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ABSTRACT 

A series of high-resolution ACDM cosniological N-body simulations are used to study the 
properties of galaxy-size dark halos as a function of global environment. We analyse halos in three 
types of environment: "cluster" (cluster halos and their surroundings) , "void" (large regions with 
density contrasts < —0.85), and "field" (halos not contained within larger halos). We find that 
halos in clusters have a median spin parameter '--^ 1.3 times lower, a minor-to-major axial ratio 

1.2 times lower (more spherical), and a less aligned internal angular momentum than halos 
in voids and the field. For masses ^ 5 x 10^^h~^MQ, halos in cluster regions are on average 
~ 30 — 40% more concentrated and have ~ 2 times higher central densities than halos in voids. 
While for halos in cluster regions the concentration parameters decrease on average with mass 
with a slope of ~ 0.1, for halos in voids these concentrations do not seem to change with mass. 
When comparing only parent halos from the samples, the differences are less pronounced but 
they are still significant. We obtain also the maximum circular velocity- and rms velocity-mass 
relations. These relations are shallower and more scattered for halos in clusters than in voids, 
and for a given circular velocity or rms velocity, the mass is smaller at z = 1 than at z = 
for all environments. At z = 1, the differences in the halo properties with environment almost 
dissapear, suggesting this that the differences were stablished mainly after z ^ 1. The halos 
in the cluster regions undergo more dramatic changes than those in the field or the voids. The 
differences in halo properties with environment are owing to (i) the dependence of halo formation 
time on global environment, and (ii) local effects as tidal stripping and the tumultuos histories 
that halos suffer in high-density regions. 

We calculate seminumerical models of disk galaxy evolution using halos with the concentra- 
tions and spin parameters found for the different environments. For a given disk mass, the galaxy 
disks have higher surface density, larger maximum circular velocity and secular bulge-to-disk ra- 
tio, lower gas fraction, and are redder as one goes from cluster to void environments. Although 
all these trends agree with observations, the latter tend to show more differences, suggesting 
this that physical ingredients not considered here as missalignment of angular momentum, halo 
triaxility, merging, ram pressure stripping, harassment, etc. play an important role for galaxy 
evolution, specially in high-density environments. 

Subject headings: cosmology: dark matter — galaxies:formation — galaxies:halos — methods:A''— body 
simulations 
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1. Introduction 

The understanding of formation and evolution 
of cosmic structures in the Universe is a fun- 
damental problem in astrophysics as well as an 
observational testbed for cosmological, particle 
and high-energy physical theories. The current 
paradigm in this topic is the inflationary moti- 
vated A- Cold Dark Matter (ACDM) hierarchi- 
cal scenario. According to this scenario, the 
backbone of galaxy formation and evolution are 
the CDM halos, which emerge from a Gaussian 
primordial density fluctuation field and assemble 
through gravitational processes. The ACDM sce- 
nario succesfuUy reproduces most of cosmological 
observational data, from the angular anisotropics 
of Cosmic Microwave Radiation to galaxy clus- 
tering (see for a review Frenk 2002). The con- 
frontation of models and observations at the level 
of galaxies is more difficult due to the complexity 
of the baryonic processes such as gas cooling, hy- 
drodynamics, star formation (hereafter SF), and 
feedback. Several approaches have been intro- 
duced to model these processes in the context of 
the ACDM scenario and their results are encour- 
aging, however, some difficulties still remain (see 
for a review Firmani & Avila-Reese 2003). To 
gain insight on galaxy formation and evolution, as 
well as on the galaxy observables that can be used 
directly to tost the ACDM scenario, one should 
identify which galaxy properties are clean tracers 
of CDM halo properties. The comparison of the 
dependences on global environment of both the 
simulated CDM halo and the observable galaxy 
properties may shed ligth on this undertaking. 

1.1. Observational dependences of galaxy 
properties on environment 

It is well known that several galaxy properties 
vary as a function of environment. The morphol- 
ogycal type of galaxies is earlier in the locally 
denser regions (morphology-density relation), be- 
ing the fraction of ellipticals maximal in cluster 
cores (Dressier 1980) and enhanced in rich (Post- 
man & Geller 1984) and poor (Tran et al. 2002; 
Tovmassian ct al. 2004) groups. An even tighter 
correlation has been found between morphology 
and cluster- or group-centric radius (Whitmore et 
al. 1993; Dommgucz ct al. 2002). The exten- 
sion of the morphology-density and morphology- 



radius relations to low local-density environment 
(cluster outskirts, low mass groups, field) has 
been a matter of debate (e.g., Postman & Geller 
1984; Domfngucz ct al. 2002). From an anal- 
ysis of SDSS data. Goto et al. (2003) have 
found that (i) in the sparsest regions both rela- 
tions flatten out, (ii) in the intermediate density 
regions (e.g., cluster outskirts) the intermediate- 
type galaxy (mostly SOs) fraction increases to- 
wards denser regions whereas the latc-type galaxy 
fraction decreases, and (iii) in the densest regions 
intermediate-type fraction decreases radically and 
early-type fraction increases. In a similar way, an 
study based on 2dF data of the luminosity func- 
tions in clusters and voids show that the popula- 
tion of faint late-type galaxies dominates in the 
latter, while, in contrast, very bright early-late 
galaxies are relatively overabundant in the former 
(Croton et al. 2005). These an other studies (e.g., 
Dommguez et al. 2001) suggest that the origin 
of the morphology-density (or morphology-radius) 
relation could be a combination of (i) initial cos- 
mological conditions and (ii) of external mecha- 
nisms (ram-pressure and tidal stripping, thermal 
evaporation of the disk gas, strangulation, galaxy 
harashment, truncated SF, etc.) that operate 
mostly in dense environments, where namely the 
relation steepens significantly. 

The morphology-environment relation evolves. 
It systematically flattens with z in the sense that 
the grow of the early-type (E-fSO) galaxy fraction 
with density becomes less rapid (Treu et al. 2003; 
Smith et al. 2004; Goto et al. 2004; Postman et 
al. 2005), the main change being in the high- 
density population fraction. Postman et al. con- 
clude that the observed flattening of the relation 
up to z ^ 1 is due mainly to a deficit of SO galax- 
ies and an excess of Sp+Irr galaxies relative to the 
local galaxy population; the E fraction-density re- 
lation does not appear to evolve over the range 
< z < 1.3. 

Observational studies show that other proper- 
ties besides morphology vary with environment 
locally (e.g., Drossier et al. 1985; Giovanelli & 
Haynes 1985; Zabludoff & Franx 1993). More 
recently it was shown that the galaxy proper- 
ties most sensitive to environment arc the inte- 
gral colour and specific SF rate (e.g., Balogh et al. 
1998, 2004; Kodamaetal. 2001; Tran et al. 2001; 
Pimbblet et al. 2002; Gomez et al. 2003; Lewis 
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et al. 2002; Kauffmann et al. 2004; Hogg ct al. 
2004; Tanaka et al. 2004), the steepest correlation 
being for intermediate mass galaxies (Kauffmann 
et al. 2004) or faint galaxies (Tanaka ct al. 2004). 
The dependences of both properties on environ- 
ment extend typically to lower local densities than 
the dependence for morphology. These properties 
are tightly related to the galaxy SF history, which 
in turn depends on internal formation/evolution 
processes related directly to initial cosmological 
conditions (e.g., Kauffmann et al. 1993; Baugh 
et al. 1996; Somerville & Primack 1999; Avila- 
Reese & Firmani 2000) as well as to external as- 
trophysical mechanisms able to inhibit or induce 
SF activity. It is difficult to disantangle the role 
of these internal or external factors, in particular 
in the densest environments. An analysis of the 
influence of environment on the CDM halo prop- 
erties and their evolution certainly would help to 
understand the role of initial conditions. 

There are also some pieces of evidence that 
the scale lenghts of spirals in clusters are sys- 
tematically smaller than those of spirals in the 
field (Aguerri et al. 2004). Besides, it was shown 
that the fraction of low surface brightness (LSB), 
blue galaxies increases toward low-density envi- 
ronments (e.g., Bothun et al. 1993; Rojas et al. 
2004), mainly in the outskirts of filamentes and in 
voids (Rosenbaum & Bomans 2004). The LSB 
galaxies may have formed in halos with high and 
well-aligned angular momentum (e.g., Dalcanton 
et al. 1997; Firmani & Avila-Reese 2000; Avila- 
Reese & Firmani 2000; Boissier et al. 2003) 
and/or have had a low SF rate history in compar- 
ison with high surface brightness galaxies due to 
their low-density environment. To explore which 
one of these processes dominate it would be help- 
ful to know if there are any differences in the halo 
angular momentum distribution and its alignment 
between voids and other environments. 

1.2. Theoretical results and expectations 

From the analysis of a 60/i~^Mpc box simula- 
tion with a particle mass of 1.57 x lO^'Mo, Avila- 
Reese et al. (1999) found that the outer density 
profiles of galaxy-size halos in cluster cores decline 
sharper and their concentrations are higher than 
for "field" halos (see also Ghigna et al. 1998; 
Okamoto & Habe 1999; Klypin et al. 1999; Bul- 
lock et al. 2001a). The last authors explored also 



the differences in concentration between distinct 
or parent halos (halos not contained within larger 
ones) and subhalos (halos contained inside larger 
halos), with the result that the latter have typ- 
ically larger concentrations than the former. As 
suggested by the referee, one should differentiate 
the effect of the "local" environment, which is sim- 
ply due to a halo being subhalo, from the global 
(large-scale) environment. Some dependence of 
the halo mass-circular velocity relation on global 
environment was also reported by Avila-Reese et 
al. (1999). 

Lemson & Kauffmann (1999) explored how 
halo mass function, formation redshift, concen- 
tration, shape and spin parameter A change with 
the density contrast of the local halo environ- 
ment. The particle mass in their simulations was 
2 X IO^^Mq or more. They found that only the 
mass function varies with environment, suggesting 
that any dependence of observable galaxy prop- 
erties on environment can be established because 
only the halo mass influences on these properties. 
Antonuccio-Delogu et al. (2002) used constrained 
initial conditions in order to produce in the same 
box two extreme environments, a void and two 
clusters (the particle mass of their simulation was 
4.2 X IO^Mq). They found some differences in the 
distribution of the spin parameter and the mass- 
velocity dispersion relation between cluster and 
void environments. 

A general conclusion of all previous numeri- 
cal works is that present-day CDM halo prop- 
erties and correlations do not strongly change 
with global environment. However, as mentioned 
above, some influence^ of cmvironment on halo 
properties and correlations has been reported. Ha- 
los in dense environments, as clusters, are ex- 
pected to collapse earlier than halos in less dense 
environments. Besides, halos in the high density 
environments suffer tidal stripping and frequent 
violent mergers. It is still an open issue the extent 
of the influence of environment on halo proper- 
tics as well as its potential effects on the baryonic 
galaxies which form within these halos. 

In this paper we will analyze high-resolution 
simulations of extreme environments, namely se- 
lected cluster and void regions, and compare 
the properties and correlations of the galaxy-size 
CDM halos from these simulations, both at red- 
shift z = and z = 1. Our aim is to explore 
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to what extent the CDM halo properties that af- 
fect observable galaxy properties change with the 
global environments. We will analyze also distinct 
galaxy-size halos ( "field" halos) and compare them 
with those in the cluster and void simulations. 

Seminumerical and semianalytical models have 
shown that most of the properties and correlations 
of galaxies formed within CDM halos depend on 
(i) the halo mass aggregation history (MAH) and 
its dominating regime (accretion or merging) as 
well as on the halo concentration that is deter- 
mined mainly by the MAH (in particular the typ- 
ical formation time) for a given mass, (ii) the halo 
spin parameter and angular momentum distribu- 
tion, and (iii) the halo mass function (e.g., Kauff- 
mann et al. 1993; Baugh et al. 1996; Mo et 
al. 1998; Avila-Reese k Firmani 2000; Firmani 
& Avila-Reese 2000; van den Bosch 2000). The 
dependence on environment of the halo mass func- 
tion has been extensively analyzed in numerical 
works (Lemson & Kauffmann 1999; Gottlober et 
al. 2003; Mo et al. 2004). Our comparative 
study here will be focused on some halo struc- 
tural and angular momentum parameters: con- 
centration, central density, shape, spin parameter, 
and the angular momentum internal alignment. 
Elsewhere we will present results related to halo 
evolution (MAH) depending on environment (see 
also Gardner 2001; Gottlober, Klypin & Kravtsov 
2001). 

In §2 we describe the method and simulations 
carried out and present the halo samples to be 
analyzed. In §§3.1 several concentration param- 
eters as well as the central density of halos from 
the different environments (samples) arc presented 
and disscused for z — and z = I. The cUiptic- 
ity distribution function of halos for the different 
samples is computed in §§3.2. Section 4 is de- 
voted to the analysis of the spin parameter and 
internal angular momentum alignment of halos as 
a function of environment. In §5 we present and 
discuss the maximum circular velocity and veloc- 
ity dispersion-mass relations of halos in the differ- 
ent environments. In §6 we discuss the implica- 
tions that the differences in the halo properties as 
a function of environment can have on the prop- 
erties of the disks formed inside these halos. A 
summary of the results and our main conclusions 
are given in §7. 



2. Method and simulations 

In our numerical simulations we adopt the flat 
cosmological model with a non-vanishing cosmo- 
logical constant (ACDM) using the following cos- 
mological parameters: Oq = 0.3, Cl\ = 0.7, and 
h = 0.7. The matter power spectrum is normal- 
ized to as = 0.9, where as is the present linear 
rms amplitude of mass fluctuations in spheres of 
radius 8/i^^Mpc. 

The simulations were run with the Adaptive 
Refinement Tree (ART) code (Kravtsov et al. 
1997). The ART code achieves high spatial resolu- 
tion by refining the base uniform grid in all high- 
density regions with an automated refinement al- 
gorithm. Initial conditions have been calculated 
either with the transfer function given by Klypin & 
Holtzman (1997) (boxes of 60 and 80 /i'^Mpc) or 
using the numerical results of the Boltzmann code 
kindly provided by W. Hu. Up to a few percent 
both transfer functions coincide and at the same 
level they coincide also with the transfer function 
provided by CMBfast. For the multiple mass sim- 
ulations we have used the mass refinement tech- 
nique described by Klypin et al. (2001). 

The halos are identified by the Bound Density 
Maxima (BDM) algorithm (Klypin & Holtzman 
1997; Klypin et al. 1999). The BDM algorithm 
first finds positions of local maxima in the density 
field. Once centers of potential halos are found, 
the algorithm identifies halos around them and re- 
moves particles which are not bound to those ha- 
los. This procedure also detects subhalos of larger 
objects - halos inside halos (for example, satellites 
of galaxies or galaxies in clusters). Particles of a 
subhalo are bound to both the subhalo and to the 
larger halo. 

Our study is focused on the analysis of the 
properties of galaxy-size CDM halos formed in 
two extreme environments: the high density en- 
vironment of clusters and the low density envi- 
ronment of voids. To this end we have at first 
performed simulations with low mass resolution. 
However, already for this simulation the initial 
displacements and velocities of the particles were 
calculated using all waves ranging from the funda- 
mental mode k = 2'it/L to the Nyquist frequency 
AiNy = 2it/LxNi/2, where A^'i = 1024 corresponds 
to the maximum possible mass resolution (cf. Ta- 
ble 1). To get the initial conditions of the low reso- 
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Table 1 
Parameters of simulations 



simulation 


Box 




Mass resolution 


Force resolution 


•^start 




[^n ivipc ) 


(tl Mq) 






Full box 












80 


80 


128^ 


2.0 X 10^" 


9.8 


50 


60A 


60 


256^ 


1.1 X 10^ 


2.0 


35 


60B 


60 


256^ 


1.1 X 10^ 


2.0 


35 


50 


50 


128^ 


5.0 X 10^ 


6.2 


60 


Multiple mass 












Clall 


80 


(5123) 


3.2 X IQS 


2.4 


50 


cir 


80 


(10243) 


4.0 X 10^ 


0.6 


50 


Ku 


80 


(10243) 


4.0 X 10^ 


1.2 


50 


^0 


50 


(10243) 


2.5 X 10^ 


0.8 


60 



lution run wc merge 83 particles to a more massive 
one. After running an inexpensive simulation with 
the resulting 1283 particles wc can select regions 
of interest. Then we use the original sample of 
small-mass particles to resimulate the regions of 
interest using the multi-mass technique described 
in Klypin et al. (2001). 

As high density regions we have selected spheres 
centered on 14 clusters. The radius of these 
spheres was typically about three virial radii of the 
cluster. Thus a significant fraction of resimulated 
halos are not within the virial radius of one of the 
clusters, but these halos are in the high-density 
environment around clusters. Voids have been se- 
lected as described in Gottlobcr ct al. (2003). 
The radius of resimulation area centered on voids 
has been taken typically about 10% smaller than 
the estimated void radius to be sure that no mas- 
sive objects around the low density region will en- 
ter the region with high mass and force resolu- 
tion, otherwise all the integration time would be 
allocated to such a massive object. We have per- 
formed high resolution simulations of selected re- 
gions in a 80 h~^Mpc box and in a 50 h~^Mpc box 
(simulations 80 and 50, respectively). For com- 
parison purposes, we have also analyzed a sample 
of galaxy-size "field" halos from two 60ft.~'^Mpc 
box simulations with the same cosmological pa- 
rameters as simulation 80 (simulations 60A and 
60B). The parameters of the simulations are sum- 
marized in Table 1. The brackets around the par- 



ticle numbers denote that all low mass particles 
were merged outside the regions of interest. 

The "field" halos in the GOhr'^Mpc box arc 
those called distinct or parent, i.e. halos which 
are not contained inside larger halos (Avila-Reese 
et al. 1999; Bullock et al. 2001a). Most of the 
parent halos are truly isolated (~ 80%); the rest 
are mostly 'binary' halos, i.e. with companions of 
masses larger than 1/3 the mass and within dis- 
tances up to 3 times the radius of the given halo. 

2.1. The halo samples according to the en- 
vironment 

From the simulations described above, we have 
five independent subsamples of halos (see Table 
2): two cluster subsamples selected from simu- 
lation 80, one made up of halos from the very 
high-resolution run Gig'', and the other built from 
the sum of the halos belonging to the 14 clusters, 
Clall ; both subsamples constitute the CLUSTER 
sample, where ~ 55% of the galaxy-size halos are 
parent (those outside the clusters), and ~ 45% 
are subhalos (mostly those within the cluster virial 
radii); two void subsamples, one consisting of ha- 
los from the simulation V^q , and the other built 
of halos from the 5 voids resimulated from the 80 
simulation, 14ii; boths subsamples constitute the 
VOID sample, where more than 90% of the galaxy- 
size halos are parent ones; one FIELD halo sample 
composed of two subsamples of galaxy-size parent 
halos from both the 60A and 60B simulations. 
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Column (2) of Table 2 gives the number of re- 
gions (clusters or voids) in each (sub)sample res- 
imulated with high resolution; this does not ap- 
ply for the FIELD sample, where all the par- 
ent galaxy-size halos in the boxes are taken into 
account. The number of halos with more than 
500 particles in each one of the samples is given 
in column (3). We set an upper limit mass of 
3 lO^^/i~^M0 for that we define as a galaxy-size 
halo. In column (4) we give an estimate of the 
density contrast [(pregion — Pu)l Pu] of our different 
environments. For the samples where several re- 
gions were used (Claii and T4n), the mean of the 
density contrasts and the standard deviation are 
reported; for the "field" halos, the density contrast 
is around zero. The mean halo mass {Np > 500 
particles) and the median spin parameter of each 
sample are given in columns (5) and (6), respec- 
tively. 

Finally, we note that in multimass simulations 
we discard the halos contaminated by massive par- 
ticles if their fraction is larger than 5% of the to- 
tal halo mass. Also fake halos (recognized by a 
too low central density or flat or decreasing den- 
sity profiles toward the center) have been removed; 
the criteria to remove these halos are a density in 
the central bin below ~ 5 lO^^h^MgMpc"^ and/or 
similar or increasing with radius densities in the 
the first 5 radial bins. To study the spin parame- 
ter distribution we demand a halo to have a mini- 
mum of 500 particles. For concentration and cen- 
tral density determinations we demand more than 
2500 particles. 

3. Structural properties 

The virialized CDM halos present a diversity of 
spherically averaged density profiles (Avila-Reese 
et al. 1998, 1999; Jing 1999; Tasitsiomi et al. 
2004) that cluster around the Navarro, Frenk & 
White (1997; hereafter NFW) fit. High-resolution 
simulations show that in the center the density 
profile is typically steeper than the NFW profile 
(Moore et al. 1999; Klypin et al. 2001; Power et 
al. 2003; Navarro et al. 2004). Unfortunately, 
there is not still a clear theoretical understanging 
of the origin of the CDM halo density profiles (see 
for some interesting results Manrique et al. 2003, 
and more references therein). Several works have 
shown that the shape of the profiles or their con- 



centration parameter if they are fitted by a NFW 
profile -which is not always a good fit, depend on 
the halo MAH (Avila-Reese et al. 1998; Firmani 
& Avila-Reese 2000; Wcchslcr et al. 2002) or even 
on the cosmological initial conditions (Avila-Reese 
et al. 2003). It was shown in Avila-Reese et al. 
(1999) that the outer density profile of CDM halos 
change with environment. Following, we will ex- 
plore whether or not the halo concentrations, cen- 
tral density, and internal shape alignment change 
with environment. 

3.1. Concentrations and inner density 

We constriict spherically averaged density pro- 
files for all the halos studied here. We can roughly 
characterize the profiles by a concentration pa- 
rameter only, because the profiles have approxi- 
mately similar shapes. The use of one parameter 
allows us to attain an easy statistical comparison 
between VOID, CLUSTER and FIELD samples. 
In the CDM hierarchical scenario one expects that 
the concentration of halos will be typically higher 
for less massive halos, because they assemble most 
of their mass earlier than larger mass halos (e.g., 
NFW). 

For the NFW profile, the concentration is de- 
fined as Cnfw = Rvir/fs, where i?vir is the virial 
radius and rs the radius where dlnp/dlnr = —2 
(NFW). The virial radius i?vir is defined as the 
radius where the average halo density is A times 
the background density according to the spheri- 
cal top-hat model. The parameter A depends on 
epoch and cosmological parameters (fioi^A); for a 
flat ACDM model, A ~ 337 and 203 at z = 
and z = 1, respectively. As mentioned above, 
not all halos can be well fitted by a NFW pro- 
file. One may introduce a concentration param- 
eter that does not depend on a particular fitting 
profile, for example the C1/5 concentration defined 
as Ci/5 = i?vir/''i/5, where ri/5 is the radius where 
1/5 of the virial mass is contained. 

For both the Cnfw and C1/5 concentrations one 
needs to measure i?vir. On one hand, this is not 
possible for a large fraction of subhalos because 
they are typically truncated at a radius i?tr < -Rvir- 
On the other hand, both cnfw and Ci/5 trace 
the global halo concentration, while galaxies form 
only in the inner parts of the halos. Thus, it is 
desirable to introduce a concentration parameter 
whose definition is independent of the virial ra- 
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Table 2 
Samples of halos 



sample 


N ■ 

^ ' regions 


nh{> 500) 


contrast 


<Mh> 
(ft-iMo) 


med 


Clall 

FIELD 


14 
1 
5 
1 


421 
318 
333 

249 
2990 


37.8 ± 18.2 

44.7 
-0.85 ± 0.05 
-0.85 
~ 


5.3 X 10" 
1.0 X 10" 
4.9 X 10^° 
2.2 X 10^° 
2.2 X 10^2 


0.024 
0.024 
0.030 
0.033 
0.033 



dius, and which measures concentration in the in- 
ner/intermediate halo regions. We introduce here 
the concentration parameter c^. To calculate c^, 
we find the radii within which the halo average 
density is equal to 4 x lO** and 4 x 10^ times the 
background density, respectively, and compute the 
enclosed mass within these radii. The ratio of 
these masses is c^. For NFW profiles in the range 
of galaxy-size halos, the typical radii where the 
halo overdensities become 4x 10^ and 4x 10^ are in 
the range of ~ 0.10— 0.45i?vir, respectively, tracing 
therefore the inner to intermediate regions. For 
example, for cnfw=12, these radii are 0.11 and 
0.35i?vir, respectively, and the ratio of the masses 
enclosed within these radii is C5=0.33. We define 
formally the truncation radius i?tr as the radius 
where the outer density profile begins to systemat- 
ically flatten out or increase {dp{r)/dt > 0). The 
radius of the halo, Rh is defined as the minimum 
beetwen i?vir and i?tr- The mass enclosed within 
Rh is the halo mass M^. 

According to our resolution tests, halos with 
more than ^ 2500 particles are needed in or- 
der to estimate reliable concentration parameters. 
This leaves us with halos more massive than ~ 
lO"/i-iM0 and - 8 X IO^/i-iMq in the void Kn 
and cluster Claii samples, respectively. Since the 
voids are defined in such a way that they contain 
only halos less massive than ^ 4 x IO^^H^^Mq 
(Gottlober et al. 2003), the ranges of masses in 
both samples do not overlap and our comparisons 
should be based on extrapolations of the mass- 
concentration dependences. However, we have one 
cluster resimulated with the same resolution as the 
voids, Clg''. Therefore, in this simulation the ha- 
los resolved with more than 2500 particles have 



masses larger than ~ lO"'^^ft.~^M0. There is also 
a void simulation, V50' , for which halos resolved 
with more than 2500 particles have masses larger 
than ~ 5 10^°h~^MQ, allowing us to expand the 
range of masses of void halos. 

3.1.1. Results at z=0 

Upper panel of Fig. 1 shows cnfw vs M^ir for 
the VOID (circles), Claii (crosses), Clg''(squeletal 
triangles) and FIELD (dots) non-truncated ha- 
los with an acceptable NFW profile fit. Non- 
truncated halos are those that attain the virial 
radius without any systematical flattening or in- 
creasing of the outer density profile. They are 
mostly parent halos. The quality of the NFW 
fit is evaluated roughly through the variance of 
the fit divided by the weighted average of the in- 
dividual variances of the data (each radial bin is 
assigned a weight according to the number of par- 
ticles inside it). Since the fit is in logarithmic 
variables, the root square of this quantitiy can be 
interpreted as a measure of the logarithmic dif- 
ference between the model fit and the data. The 
latter can be expressed as a percentage deviation; 
from visual inspections, we have found that an ac- 
ceptable NFW fit to the density profile is when this 
deviation is less than ~ 15%. More than 80% of 
the non-truncated halos obey this condition. The 
parent halos in the CLUSTER sample are high- 
lighted with an open square. As expected, most 
of the non-truncated halos are parent ones. In the 
VOID sample, more than 95% of the halos with 
more than 2500 particles are parent ones, then, 
we do not highlight these halos in the figures in 
order to avoid overplotting. The FIELD sample 
has been constructed only with parent halos. 



7 



In the lower panel of Fig. 1, we show Ci/5 vs. 
A/vir using the same symbols than in the upper 
panel. We recall that the definition of Ci/5 is not 
related to any postulated density profile. There- 
fore, in the lower panel halos with non-acceptable 
NFW fits are also included. In Fig. 2 we present 
cg vs Mh (upper panel), and p_2 vs Mh (lower 
panel) with the same symbol code as in Fig. 1. 
The parameter p_2 is a characteristic inner den- 
sity measured at the radius where the slope of the 
logarithmic density profiles becomes —2. The def- 
initions of the concentration cg and of p_2 do not 
depend on i?vir, therefore, in Fig. 2 both trun- 
cated and non-truncated halos are included. 

The concentration parameters and the central 
density systematically decrease with mass for the 
CLUSTER sample. In all the cases the scatter is 
large and the correlations are weak. The Pearson 
correlation coefficients of the linear regressions (in 
logarithmic variables) are between -0.35 and -0.55. 
As a reference, in Table 3 the parameters of the 
direct linear regressions are given. We also report 
the global dispersion of each linear regression, a. 
If the scatter is Gaussian distributed (lognormal 
in linear variables) and independent of mass, then 
a represents the la (68%) variation of A(logy) at 
any mass, where Y is one of the concentration pa- 
rameters or p_2- By multiplying cr by 2.3, one get 
the width of the lognormal distribution (in natural 
logarithm), ain ~ 2.3a. 

Solid lines in Figs. 1 and 2 indicate the lin- 
ear regresions to the Claii+Cl[?''(CLUSTER) sam- 
ple data. The slopes of the regressions for cnfw 
and Ci/5 are shallow. If we take into account all 
the halos, and not only those that are well fitted 
by the NFW profile, then the slope for logCNFW 
vs logMvir is steeper, —0.14, but the scatter be- 
comes larger (from cr = 0.13 to 0.17). Note that a 
significant fraction of the halos in the cluster sim- 
ulations (^ 55%) are outside the virial radii of the 
cluster halos, i.e. they are not subhalos, but are 
in any case in overdense regions. 

The cs parameter is independent of the defini- 
tion of the halo radius (virial or truncation); in this 
sense, cs is a more direct estimator of halo mean 
(inner) concentration, and it allows to include al- 
most all the subhalos in the sample. As Fig. 2 and 
Table 3 show, for halos in CLUSTER sample, cs 
on average decreases as Mh increases. This con- 
firms the decreasing of the physical concentration 
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Fig. 1. — The NFW and 1/5 concentration param- 
eters vs. virial mass A/virfor the Claii, Clg" , VOID 
(K1II+ Vio ), and FIELD samples. The symbols are 
explained in the lower panel. Halos with truncated 
radius are not included in the plots. In the up- 
per panel were included only halos whose density 
profiles are fitted well by the NFW model. Solid 
lines are direct linear regressions to the CLUSTER 
(Claii-f Clg'') sample, while the error bars are the 
dispersions of the VOID sample (the abscise is ar- 
bitrary since there is no dependence on mass for 
this sample). Thick dashed line in the upper panel 
is the prediction of a toy model calibrated to N- 
body simulations of relaxed isolated halos (Eke, 
Navarro, & Steinmetz 2001). 

of halos/subhalos in dense environments as their 
masses (virial or truncated) increase. 

The dispersions of the linear regressions pre- 
sented in Table 3 give an estimation of the scatter 
of the concentration parameters and p-2. As Bul- 
lock et al. (2001a) found, the scatter for Cnfw 
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Fig. 2. — Inner concentration parameter, c^, and 
P-2 density vs. halo mass Af/,. for the Claii, Clg" , 
VOID (Van+ ViS'), and FIELD samples. The sym- 
bols are explained in the upper panel. Solid lines 
are direct linear regressions to the CLUSTER sam- 
ple, while the error bars are the dispersions of the 
VOID sample (the abscise is arbitrary since there 
is no dependences on mass for this sample). Thick 
dashed line in the lower panel is the prediction of 
a toy model for the NFW profile calibrated to N- 
body simulations of relaxed isolated halos (Eke, 
Navarro, & Steinmetz 2001). 

does not depend significantly on mass and it has 
a lognormal distribution. These authors report 
that the la deviation of A(logCNFw) is 0.18 for 
distinct halos, and 0.24 for subhalos (the correct 
values are actually 0.14 and 0.18, respectively, J. 
Bullock, private communication). In the under- 
standing that the dispersion cr reported in Table 3 
is similar to the Icr deviation of A(logCNFw) (see 
above), we find that the scatter of Cnfw for our 



prunned CLUSTER sample (including both par- 
ent halos and subhalos) is similar to that of the 
distinct halos at all environments in Bullock et al. 
(2001a); by including halos with density profiles 
not well fitted by the NFW model (as in Bullock 
et al. 2001a), the scatter of the CLUSTER sam- 
ple increases to 0.17. By separating the prunned 
sample in parent halos and subhalos, the la vari- 
ations of A(logCNFw) are 0.12 and 0.17, respec- 
tively. We have calculated also the value of the 
1—a deviation of A(logCNFw) in a range of masses 
from 5 10" - 2 IO^^Mq for the halo/subhalo sam- 
ple; the value is the same than that found through 
the linear regression to all the mass range, 0.13. 

For the VOID (low mass) halo sample we do 
not find any trend of cnfw, C1/5, c^, and p_2 with 
mass, although the range of masses studied is too 
small to claim for definitive conclusions. It seems 
that the halo assembling process in low-density en- 
vironments is so slow, constant and late that a 
posible mass-formation epoch relation (Lacey & 
Cole 1993; Avila-Reese et al. 1998) loses sense; 
the dependence of concentration on mass is at- 
tributed to this relation (e.g., NFW). In the case 
of halos in VOID sample, almost all of them are 
parent ones. The error bars in the left sides of 
Figs. 1 and 2 represent the (logarithmic) la dis- 
persions of Cnfw J Ci/5j Ci5 and p_2 for the VOID 
sample (see also Table 3 for roughly the same re- 
sults) . Since there is no dependence of these quan- 
tities on mass, the abscise of the bars is arbitrary. 
The scatters of the concentrations and /0_2 in the 
VOID halo sample are significantly smaller than 
the ones in the CLUSTER sample. For example, 
the la variation of A(logCNFw) is 0.09 for the for- 
mer sample compared to 0.13 for the latter one. If 
one takes into account only the parent halos, then 
the scatters are 0.09 and 0.12, respectively. There- 
fore the difference in the scatter of Cnfw between 
the VOID and CLUSTER halo samples is mainly 
due to a (global) environmental effect. The same 
applies for the scatters of the other concentration 
parameters and p-2- 

The concentrations and p-2 of the FIELD (par- 
ent) halos tend to decrease with mass, but the 
trends are shallower than for the CLUSTER halo 
sample. We notice that the trends of logCNFW and 
logCi/5 with logAfvir for our FIELD sample are 
shallower than those reported for isolated halos in 
previous works, where erg = 1 was used (Avila- 
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Reese et al. 1999; Bullock ct al. 2001a), at 
least in the range of masses explored here. In 
Figs. 1 and 2 we show a toy-model prediction 
for cnfw and p_2 (long-dashed lines) following 
Eke, Navarro, & Steinmetz (2001) for the same 
cosmological parameters than in our simulations 
and (Tg = 0.9. The toy model for cnfw was nor- 
malized to a set of isolated halos in a wide range 
of masses (re)simulated with very high resolution 
(Eke, Navarro, & Steinmetz 2001). The toy model 
roughly agrees with our FIELD halo results. 

Interestingly, studies aiming at infering the halo 
Cnfw concentration from galaxy observables show 
that the slope of the Cnfw — -^vir relation is signif- 
icantly steeper for halos of elliptical galaxies (e.g., 
Sato et al. 2000), which arc localized typically 
deep inside the clusters, than for halos of spiral 
galaxies (Jimenez, Verde & Oh 2003), which are 
localized in less dense environments. 

Our main interest here is to compare the prop- 
erties of halos in low and high-density environ- 
ments. From Figs. 1 and 2, one concludes that 
the concentrations and p^2 of the VOID sample 
are lower on average than those of the CLUS- 
TER sample, although the scatter in both sam- 
ples is large. In fact this conclusion is based 
mostly on extrapolations of trends because the 
VOID and CLUSTER samples almost do not over- 
lap in mass. At ^ 10^^/i~^Mq, the concentrations 
Cnfw J and c^ are on average 40% smaller for 
halos in the VOID sample than for halos in the 
CLUSTER sample, and the inner density p^2 is 
on average two times lower for the former than for 
the latter. For larger masses, all these differences 
reduce. There arc also differences in the scatters: 
they are larger for the CLUSTER sample than for 
the VOID and FIELD ones (see Table 3). 

The halos from the FIELD sample also tend 
to have on average lower concentrations and cen- 
tral density than those halos from the CLUSTER 
sample, although the differences arc smaller than 
those between VOID and CLUSTER halos, in par- 
ticular for /9_2. If we assume that the prediction 
for Cnfw given by the Eke, Navarro, & Steinmetz 
(2001) toy-model is a good fit to our FIELD sam- 
ple data, then its extrapolation to lower masses 
shows that the difference in cnfw between FIELD 
and CLUSTER halos remains roughly the same 
(Fig. 1). For p_2, the difference seems to increase 
at lower masses (Fig. 2). 



have also compared only the parent halos 
from the CLUSTER sample with those from the 
VOID and FIELD samples. The differences in 
concentrations and p_2 become less pronounced 
than in the case when the subhalos were included. 
The slopes of the cnfW) Ci/5, C5 and p-2 vs mass 
relations for the CLUSTER parent halos reduce 
to -0.07, -0.06, -0.07 and -0.16, respectively (to 
be compared with those given in Table 3), be- 
coming closer to those of the FIELD sample, but 
still steeper. As expected, c^ changes more than 
Cnfw and ci/5 because the definition of the lat- 
ter concentrations imply that the halos should be 
non-truncated, and the non-truncated halos are 
mostly parent ones, so that the fraction of CLUS- 
TER subhalos in Fig. 1 is small in any case. We 
have also compared the concentrations and p-2 
of the very few subhalos from the VOID sam- 
ple with the extrapolations of these parameters 
to the corresponding masses from the CLUSTER 
subhalo sample. All the VOID subhalo concentra- 
tions and p_2 he below the corresponding extrap- 
olations from the CLUSTER subhalos. 

Our results show that halos in dense environ- 
ments are on average more concentrated, with 
higher central densities and with larger scatters 
in these parameters than halos in low-density en- 
vironments. We find that these differences are ow- 
ing to both nature and nurture reasons. Related 
to the former is mainly the fact that halos in dense 
regions typically assemble most of their masses 
earlier than halos in low dense regions, through 
violent MAHs. Sheth & Tormen (2004) indeed 
have shown that halos in dense regions form at 
earlier times than do halos of the same mass in 
less dense regions. Related to "nurture" is the ef- 
fect of steepening of the outer (sub) halo profile 
due to tidal stripping by the parent halo (Ghigna 
et al. 1998; Okamoto & Habe 1999; Avila-Reese 
et al. 1999: Klypin et al. 1999). This effect 
tends to make the subhalos more concentrated and 
it is subject to a local condition rather than to 
the global environment. Subhalos are indeed sys- 
tematically more concentrated than parent halos, 
however, both increase also their concentrations as 
the environment density increases (see also Bul- 
lock et al. 2001a). The effect of the formation 
epoch, related to the global environment, affects 
likely more the central halo density than concen- 
trations, while the inverse is expected for the local 
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evolutionary effects. As will be seen below, the 
differences with environment from 2; = 1 to z = 
are more pronounced for concentrations than for 

P-2- 

3.1.2. Results at z^l 

We analyze the CLUSTER, VOID and FIELD 
samples at z = 1 in the same statistical sense as 
done at z = 0. Similar to Figs. 1 and 2, we show 
in Figs. 3 and 4 the concentrations and p-2 versus 
mass but now at z = 1, and in Table 3 the param- 
eters of the corresponding linear regressions are 
given. The first impresion is that the differences 
in the concentrations and central density seen at 
z = among the three halo samples are hardly 
present at z = 1. Our results thus show that at 
this epoch the Cnfw, C1/5, c^ and p-2 versus mass 
relations do not depend on environment. 

The central concentration p_2 for the CLUS- 
TER halos is on average only a little larger at 
z = 1 than at z = 0. Halos in dense environ- 
ments assemble their present-day mass early in 
such a way that their central densities are stab- 
lished likely before z = 1. For the VOID halos, 
the mean of p-2 at z = 1 is larger than the one 
at z = by ~ 70%. Halos in low-density environ- 
ments assemble slowly, incorporating most of their 
mass lately. 

The concentration parameters change more 
with z than the central density. As previously re- 
ported (e.g., Bullock et al. 2001a; Eke, Navarro, 
& Steinmetz 2001; Navarro et al. 2004), the 
Cnfw concentration parameter is lower at higher 
redshifts for a given mass. In the upper panel of 
Fig. 3, the dashed curve shows the cnfw ^ A/vir 
dependence ai z — 1 given by the toy model of 
Eke, Navarro, & Steinmetz (2001). According 
to this model Cnfw (1 + z)^^ approximately 
(see also Bullock et al. 2001a). In Figs. 3 and 
4 are also shown (thin solid lines) the linear re- 
gression for the CLUSTER sample and the la 
scatter of the VOID sample plotted in Figs. 1 and 
2 (z = 0). The concentration cnfw is on aver- 
age ^ 1.5 and 2.1 times lower at z = 1 than at 
z = for the VOID and CLUSTER samples, re- 
spectively. For the latter sample, the comparison 
is made at IO^^/i^^Mq. For the FIELD sample, 
the differences of Cnfw between z = and z = 1 
are on average slightly less than a factor of two, at 
least in the range of masses explored here. Thus, 
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Fig. 3. — Same as in Fig. 1 but for a snapshot 
at z = 1. Thick lines and error bar correspond to 
the current samples (z = 1), while thin line and 
error bar are the same ones as shown in Fig. 1 
(z — 0). The bars were slightly shifted to avoid 
overlapping. 

the evolution of Cnfw seems to be slightly dif- 
ferent in the different environments. A similar 
behaviour is seen for Ci/5 and c^. The scatters in 
all the measured parameters are larger at z = 1 
than at z = 0, reflecting likely the fact that at 
z = 1 the halos are in general less relaxed than at 
z = 0. 

The increase of concentration with time is more 
pronounced in the cluster environment. We ana- 
lyze the behaviour of concentration for halos inside 
the present-day cluster virial radius (subhalos) in 
simulation Clg''. Because for this simulation we 
have the halo MAHs, concentrations for a given 
halo can be measured at z = and z = 1. We 
find that the ratios of concentrations measured 
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Table 3 

Linear regression parameters (concentrations and p_2) 



Sample IgCNFW = a + WgAfvir lgci/5 = c + rflgMvir lgC5 = e + /lgMft lgp-2 = g + hlgMh 
aba c d a e f a g h a 



z = 


























CLUSTER 


2.23 


-0.09 


0.13 


2.10 


-0.09 


0.11 


0.77 


-0.11 


0.11 


19.01 


-0.29 


0.33 


VOID 


1.10 


0.00 


0.09 


0.96 


-0.00 


0.07 


-0.38 


-0.01 


0.07 


15.04 


0.04 


0.29 


FIELD 
z = 1 
CLUSTER 


1.51 


-0.04 


0.11 


1.53 


-0.05 


0.08 


0.00 


-0.05 


0.10 


15.94 


-0.05 


0.32 


1.40 


-0.05 


0.15 


1.36 


-0.05 


0.09 


-0.23 


-0.04 


0.18 


17.47 


-0.15 


0.40 


VOID 


1.15 


-0.02 


0.10 


1.22 


-0.04 


0.06 


-0.34 


-0.02 


0.13 


17.30 


-0.14 


0.40 


FIELD 


0.67 


0.00 


0.15 


0.60 


0.01 


0.09 


0.15 


- 0.04 


0.19 


17.17 


-0.14 


0.38 



at 2: = and z = 1 tend to increase as the 
Mh{z = Q)/Mh{z = 1) ratio decreases (Fig. 5); 
i.e., halos that have suffered more mass loss due 
to tidal stripping, end preferentially with higher 
concentrations. Therefore, as mentioned above, 
tidal stripping seem to be an efficient mechanism 
for increasing the average concentration of halos. 
Moreover, since the concentration is an increasing 
function of the mass loss, then some steepening 
of the concentration-mass relation in high-density 
environments is expected. This is what we ob- 
serve in our simulations, in particular for cs and 
Ci/5 (Figs. 3 and 4). For cnfw this effect is not 
seen because only halos that have a good fit to the 
NEW profile are selected, and these are typically 
well relaxed halos that have not felt strong tidal 
stripping. 

3.2. Halo shape distribution 

Halo axis ratios as well as the directions of the 
principal axes are determined by iteratively^ diag- 
onalizing the tensor 



lij = 51 a;^a;J•/a^ a= ^x^ + ^ + —. (1) 

The sum is over all particles within the halo ra- 
dius Rh (virial or truncation radius, whichever 

'The first iteration consider all particles inside a sphere of 
radius equal to the halo radius. For the second iteration, 
we reorient the system so as the x axis lies on the major 
axis of the ellipsoid found in the first iteration. Wc now 
use all particles inside the ellipsoid. The iteration proceed 
untill convergence is obtained. 



is smaller), xt {i = 1, 2, and 3) are the parti- 
cle coordinates with respect to the halo center of 
mass, s is the short- to-long axis ratio, and q the 
intermediate-to-long axis ratio. We use here halos 
with more than 500 particles. 

To obtain an estimate of the error in the de- 
termination of the axial ratios due to low particle 
number, we analyze the well-resolved chistcir halos 
from our simulations. For each one of these halos, 
which contain hundreds of thousands of particles, 
we randomly extract only 500 particles and mea- 
sure q and s in 100 realizations. We found the the 
la errors in q and s amount to about 10% and that 
they do not depend on the mass of the cluster-size 
halo. These errors are an upper limit because our 
galaxy-size halos have more than 500 particles. 

Figure 6 shows the cUipticity as a function of 
halo mass at two redshifts for the FIELD halo 
sample. For a given epoch, low-mass halos are 
preferentially more spherical than high-mass ha- 
los while, for a given mass, halos at z = 1 are on 
average more flattened than halos at z = 0. These 
results agree with those of Bullock (2001) and 
Hopkins et al. (2005) and seem to suggest that 
present-day halos, having had more time to relax 
and attain a more isotropic orbit distribution, as 
compared with high-redshift ones, tend to be more 
spherical. Moreover, since small halos form earlier 
on average than big ones, they also have had more 
time to relax and thus become more spherical. 

In Fig. 7 we show the distributions of the halo 
ellipticity ei = 1 — s for our different samples. 
The ellipticity distribution of sample Claii is shown 
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Fig. 4. — Same as in Fig. 2 but for a snapshot at 
z — 1. According to the spherical coUapse model, 
the overdensity is 1.65 times smaller at z = 1 than 
at z = for our cosmology; this was taken into ac- 
count in cs- Thick lines and error bar correspond 
to the current samples {z = 1), while thin line and 
error bar are the same ones as shown in Fig. 2 
(z = 0). 

in left panel while the one for sample VOID is 
shown in right panel. They are both well fit by a 
lognormal distribution, 



Pie) 



1 



ev27rcr 



exp 



2a2 



(2) 



with best fit values (cr,eo) given on the correspond- 
ing panels. Curves on panels are the lognormal 
best fits to the data. The lognormal fit of sample 
FIELD is shown as a dotted line in both panels 
for comparison purposes. 

A Kolmogorov-Smirnov (KS) test shows that 
the probability that Claii and FIELD halo sam- 
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Fig. 5. — Ratio of the concentration parameter 
Cs measured at z = and z = 1 vs the ratio 
of the halo masses at these epochs for individ- 
ual halos inside the present-day virial radius of 
cluster Clg''. Halos with more than 500 parti- 
cles at z = are shown. Those with more than 
2500 particles are marked with solid circles. Ha- 
los that have lost more mass, i.e. that underwent 
more tidal stripping, have systematically increased 
more their concentrations than halos less affected 
by tidal stripping. 

pies are drawn from the same parent distribution 
is very small, 3.0 x 10^^. Thus, according to 
Fig. 7, halos in clusters and their outskirts are 
on average more spherical than halos in voids or 
in the field. This conclusion may result at first at 
odd because, being halos in clusters more massive 
on average than their counterparts in voids, they 
are expected to be more fiattened (see Fig. 6). 
However, the environmental effect related to the 
halo formation epoch seems to dominate over the 
weak dependence of ellipticity on mass. Halos in 
a high-density environment assemble earlier than 
in a low-density one, and as follows from Fig. 6, 
halos tend to be more spherical as they assemble 
earlier. On the other hand, when the comparison 
is made between halos from the VOID and FIELD 
samples, the mass effect seems to dominate over 
the environmental one, in such a way that the for- 
mer are on average less flattened than the latter. 
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Fig. 6. — The average ellipticity, e = 1 — s, as 
a function of halo mass at two redshifts z = 
and 1 for halos in the FIELD sample. Error bars 
reflect the Poisson uncertainty associated with the 
number of halos within the bin and not the scatter 
in the relation. 
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Fig. 7. — Ellipticity distribution (histograms) for 
halos with more than 500 particles drawn from a 
cluster environment (left panel) and from a void 
environment (right panel). Curves on each panel 
are lognormal best-fits to the data. The lognor- 
mal fit of the e distribution for the FIELD sample 
of halos is shown in both panels for comparison 
purposes. The values of the parameters of the fits 
(<Te, eo) are also shown in panels. 

4. Spin parameter and angular momen- 
tum alignment 

4.1. Spin parameter distribution 

We compute the halo total angular momentum 

as 

n 

J = ^m^ri X Vi, (3) 



where and Vi are the position and velocity of 
the ith particle with respect to the halo center of 
mass. We follow Bullock et al. (2001b) and define 
a modified spin parameter A' to characterize the 
global angular momentum of a halo 



V2MhVhRh 

where Jh is the angular momentum inside the halo 
radius Rh (see §3.1), and Vh is the circular veloc- 
ity at Rfi- Hereafter, in this section, we drop the 
prime in A'. 
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Fig. 8. — Spin parameter distribution (his- 
tograms) for halos with more than 500 particles 
drawn from a cluster environment (Claii sample, 
left panel) and from a void environment (VOID 
sample, right panel). Curves on each panel are 
lognormal best fits to the data. The lognormal fit 
of the A distribution for halos in the FIELD sample 
is shown in both panels for comparison purposes. 
The values of the parameters of the fits (Ao,o'a) 
are also shown in panels. 

Figure 8 shows the spin parameter distributions 
for our different galaxy-size halo samples. In this 
case, halos with more than 500 particles were used. 
The halos that are subhalos are cut at the trun- 
cation radius. The spin distributions of the Claii 
and VOID samples are shown in left and right pan- 
els, respectively. Curves on panels are the corre- 
sponding lognormal best fits to the data. The log- 
normal fit of sample FIELD is shown as a dotted 
line in both panels for comparison purposes. The 
(ct,Ao) parameters of the lognormal distributions 
are given inside the panels. The median values of 
A for the different samples are presented in Table 
2. The conclusion from Fig. 8 and Table 2 is clear: 
halos from the Claii sample have on average lower 
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spin parameters and a wider distribution than ha- 
los from the VOID and FIELD samples. The Ao 
and fT parameters of the former sample are 30-40% 
and ~ 15% times smaller than for the latter sam- 
ples, respectively. 

The difference between the A distributions of 
halos in the Claii and VOID samples is confirmed 
by a KS test: the probability that both sam- 
ples belong to the same parent distribution is 
1.17 X 10^'"'. Our result agrees with that obtained 
by Reed et al. (2005). They also found that halos 
in high-density environments have smaller spins 
than halos in the field. Arc these differences due 
to global environment effects or to the fact that 
CLUSTER subsamples have a higher fraction of 
subhalos? To address this question we compare 
the A distributions of the parent halos and sub- 
halos of the CLUSTER subsamples. By applying 
a KS test we find only slight differences in these 
distributions: the probabilities that both samples 
belong to the same parent distribution are 0.89 for 
the Clg'' sample and 0.45 for the Claii sample. On 
the other hand, the A medians (or Aq) of the halo 
and subhalo samples in both cases agree within 
5%. We notice that the FIELD sample is by con- 
struction composed of only parent halos, and the 
VOID halos are esentially parent ones (more than 
90%). 

How much of the difference we have found can 
be attributed to cluster-to-cluster scatter? Unfor- 
tunately, we only have one cluster, Clg"^, with hun- 
dreds of halos, each one with more than 500 parti- 
cles inside its virial or tidal radius. The other clus- 
ters are less resolved and they thus have less halos 
than Clf" by far. We avoid the halo low-numher 
problem by combining these clusters to form only 
one, which in principle is independent from Clg". 
We measure P{X) for the independent Clg"^ halo 
sample and test whether this distribution is sim- 
ilar to the one obtained for the composite Claii 
sample. We found differences in the values of the 
fitting parameters (cr,Ao) of less than 5%. Slightly 
higher differences are measured when the compar- 
ison is made between the P{\) distribution from 
the same sample Claii but at two close snapshots. 

Our results show that present-day halos in 
dense environments have a spin parameter dis- 
tribution shifted to lower values with respect to 
those of halos in the field or in the voids. Does 
A decrease due to the lost of high angular mo- 



mentum material by halos that suffer strong tidal 
stripping in dense regions? A way of testing this 
hypothesis is by simply measuring P(A) in the 
cluster environment at much earlier time, when 
most halos have yet to experiment strong tidal 
stripping. We find a median value \med = 0.036 
at z = 1 for the CLUSTER sample, which is com- 
parable to Xmed = 0.033 for the FIELD sample 
at present time. This result is consistent with the 
hypothesis above but does not actually prove it. 
On the other hand, we showed above that sub- 
halos and parent halos have, according to a KS 
test, similar P(A) distributions. These results ap- 
pear to contradict each other. Notice, however, 
that the histories of parent halos in high-density 
regions can be very different from their counter- 
parts in low-density environments because the for- 
mer undergo a "tumultuous life" . They may have 
suffered, for instance, tidal stripping in the past, 
from close encounters with major substructures, 
or their mass accretion could have been stopped 
(Kravtsov et al. 2004b). We also looked for a 
cluster-centric radial A dependence in Claii and 
Clf samples. Unlike Reed et al. (2005), we did 
not find any systematical decreasing of Xmed as 
the cluster-centric radius is smaller. Notice that 
Reed et al. (2005) analyzed subhalos with more 
than 144 particles, while our subhalos have more 
than 500 particles; many subhalos are below this 
limit, but the measure of A for them is not reliable. 

For the VOID and FIELD samples, Xmed = 
0.034 and 0.036 at z = 1, respectively. There- 
fore, A for these halo samples also decreases on 
average from z = 1 to z = 0, but very little. In 
this case, the small decreasing of A with time could 
be explained by the accretion mode proposed by 
Peirani et al. (2004). In summary, while halos in 
low-dense regions seem to keep constant or slightly 
decrease their A values from 2; = 1 to 2; = 0, ha- 
los from the cluster-like regions tend to decrease 
significantly their A values. We interpret the lat- 
ter as a consequence of two effects: (i) the earlier 
halo assembly epochs typical of higher density re- 
gions, and (ii) the subhalo tidal stripping and the 
the tumultuous mass assembly history of halos in 
clusters and their surroundings. For a more de- 
tailed discussion of these effects it is necessary to 
follow the individual mass and angular momen- 
tum assembly histories of the halos. In a subse- 
quent paper we will construct and analyze these 
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individual histories for each one of our samples. 

4.2. Internal angular momentum align- 
ment 

In previous subsection we explored the distribu- 
tion of the magnitude of the angular momentum 
as a function of environment. Now we will explore 
possible differences in the distribution of the an- 
gular momentum internal alignment. We measure 
the typical alignment by the angle 9i/2 between 
the mean angular momentum of the particles in 
the inner sphere of half-mass radius and the mean 
angular momentum of the particles in the outer 
half-mass spherical shell (Bullock et al. 2001b). 
As discussed by Bullock et al. (2001b) large er- 
rors are involved in the determination of the direc- 
tion of the angular momentum vectors and thus 
in cos 01/2. The error in cos 611/2 is estimated in 
two manners. In the first one, we assign errors 
to each component of the angular momentum, in 
both inner and outer half-mass regions, using a 
Monte Carlo procedure similar to the one built to 
estimate the error in the axial ratios (§3.2), and 
perform a standard propagation of errors. In the 
second one, the Monte Carlo procedure is applied 
directly on cos0i/2. The first method gives larger 
errors. This is so likely because the internal and 
external angular momentum are not entirely in- 
dependent quantities. We found that the error 
depends not only on the number of particles, Np, 
and A but also on the intrinsic alignment: a more 
aligned halo have a smaller error. For Np as low 
as ^ 500, the \a error in cos 6*1/2 can reach 50% of 
the measurement. Bearing this in mind, we now 
proceed to present our results (still uncertain) and 
their possible intepretations. 

The distributions of cos0i/2 for our different 
environments at 2; = are shown in Fig. 9. As 
it was the case for -P(A), we use halos with more 
than 500 particles. Figure 9 shows that the Claii 
halos are on average less aligned than the halos 
from the FIELD and VOID samples. The median 
values of cos 6*1/2 for the Claii FIELD, and VOID 
samples are 0.51, 0.80, and 0.83, respectively. In 
order to find a possible explanation for these dif- 
ferences with environment, we have also measured 
the alignment distributions at a much earlier time, 
z = 1 (Fig. 10). Figures 9 and 10 show that the 
population of VOID and FIELD halos at z = is 
more aligned on average than the one at z = 1. 



A more quantitative estimate of this difference is 
confirmed by the KS test: for the FIELD sample, 
for example, the probability that the populations 
at 2; = and z — 1 belong to the same parent 
distribution is 4.5 x 10~^, with cos0i/2^mod — 0-80 
and 0.75, respectively. Regarding the Claii sample, 
there is no any significant difference in the align- 
ment distributions of halos at z = and z ~ \. 
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Fig. 9. — Angular momentum alignment distribu- 
tions for our four samples of halos at z = 0. Two 
different distributions are shown in each panel for 
comparison purposes. Distributions shown as dot- 
ted histograms are also shaded for more clarity. 
Halos in voids and those in the field present similar 
distributions but they both are on average more 
aligned than their counterpart in clusters. 

We remark that the results presented above 
need to be confirmed by future analysis with bet- 
ter resolved halos. If confirmed, we find the 
following explanation for the evolutionary effect 
seen in the internal angular momentum align- 
ment of FIELD and VOID sample halos. On one 
hand, most of the angular momentum of the outer 
sphere, which is where essentially most of the halo 
angular momentum resides, was acquired during 
the linear grow. On the other hand, it is expected 
that most of the angular momentum of the in- 
ner sphere comes from the violent, initial merger- 
growth phase (Vitvitska et al. 2002; Peirani et al. 
2004). The material accreted by the halo at later 
times brings angular momentum oriented in the 
direction of the angular momentum of the outer 
shell. As part of this material is incorporated into 
the inner sphere the alignment increases. This 
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mechanism is not expected to apply to halos in 
clusters or their outskirts because they typically 
do not incorporate material as soon as they fall 
into the cluster halo; in some cases the halos in 
clusters even loose mass due to tidal stripping. 

Regarding the question of why present-day ha- 
los in clusters are less aligned than those in the 
field or in voids, a comparison between Claii and 
FIELD halos at z = 1 shows that they are equally 
well aligned. This is probably so because at z = 
1 these environments are still not too disimilar. 
However, unlike halos in the field and voids, ha- 
los in clusters evolve under the influence of strong 
tidal fields. Halos in sample CLUSTER have a 
mass assembly history very different from their 
counterparts in samples FIELD or VOID. In the 
cluster environment, a significant fraction of ha- 
los present at z = 1 do not survive until z = 0, 
others end up with their masses significantly re- 
duced, while some others, those that at present 
day are at cluster outskirts, may grow as what is 
typical of a halo in field and void environments. 
Figures 9 and 10 show that halos in sample Cl|'' 
at z = are more disaligned than halos ai z — 1. 
For the subset of halos of sample Clg'' that are 
within i?vir at z = (subhalos), which were also 
identified at z: = 1, we find the same trend: the 
halos at z = 1 are on average more aligned than 
the halos at z = 0. A halo by halo comparison, 
however, shows that the trend is not systematic; 
i.e., there are halos that exhibit a higher alignment 
at z = 0. In summary, it seems that in an envi- 
ronment in which halos stop growing, an increase 
in alignment is not expected, but even more, it 
seems that a cluster-like environment acts on the 
halo angular momentum internal alignment in a 
non-trivial way: the change in alignment varies 
halo by halo. 

4.3. Alignment between the halo shape 
and angular momentum 

Halo axis ratios as well as the directions of the 
principal axes are determined as explained in §3.2. 
We denote with 9 the angle between the angu- 
lar momentum axis within and the direction 
of the minor principal axis. In Fig. 11 we show 
the distributions of 6* at z = for our different 
halo samples. We use here also halos with more 
than 500 particles, although the iterative proce- 
dure sometimes need more than that to converge. 




Fig. 10. — Angular momentum alignment distri- 
butions for our four samples of halos at z = and 
z = 1. Each panel correspond to a halo sample. 
Distributions shown as dotted histograms are also 
shaded for more clarity. Halos in voids and those 
in the field are on average less aligned at z = 1 
than at z = 0. 

The number of halos used to draw the histograms 
in Fig. 11 can thus be lower than those shown in 
Table 2 (see §2). We repeat here the procedure we 
follow to evaluate the la error in axial ratios and 
cos 01/2 to obtain the corresponding la error in 
cos 9. Errors in cos 6* are smaller than those ones 
in cos 01/2- They amount to about 20% for halos 
with 500 particles. Taking into account that this 
particle number is a lower limit, we expect errors 
in cos^ to be on average less than 20%. On the 
other hand, errors in cos 9 are also dependent on 
the degree of alignment but they seem to be less 
sensitive to it than the errors in cos0i/2- 

The angular momentum axis of halos tends to 
align to the minor axis of the halo ellipsoids (in 
a plot like Fig. 11, a random oriented distri- 
bution would be uniform). This appears to be 
a generic prediction of the CDM cosmology (see 
also Faltenbacher et al. 2002; Bailin & Steinmetz 
2004). Besides, we find some trend of decreasing 
the alignment angle 9 from Claii to FIELD and 
VOID halo samples. This visual result is mildly 
supported by a KS test: the probability that the 
Claii and VOID samples belong to the same parent 
distribution is 2.4 x 10~^. 

The relative orientations of the halo angular 
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momentum direction and the halo principal axes 
could have important consequences on disk galaxy 
formation, the dynamics of the disks inside the 
halos as well as on observational studies aimed to 
explore the distribution of satellite galaxies with 
respect to the parent halo shape and orientation. 
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Fig. 11. — Distributions of the alignment of the 
minor principal axis and the angular momen- 
tum within Rh for our four samples of halos at 
z = 0. Two different distributions arc shown in 
each panel for comparison purposes. Distribu- 
tions shown as dotted histograms are also shaded 
for more clarity. Halos in voids and those in the 
field present similar distributions but they both 
are on average more aligned than their counter- 
part in clusters. 



5. Structural relations 

5.1. Maximum circular velocity vs mass 

Figure 12(a) shows maximum circular veloc- 
ity, Knax, vs. mass, M/i, for the halos from all 
the samples analyzed here. As usual, Knax is 
the maximum value of the smoothed radial pro- 
file Vciic('') = y/GM{r)/r, computed under the 
assumption of spherical symmetry. The solid and 
dashed lines are the linear regressions for the 
CLUSTER and VOID samples, respectively. The 
Pearson correlation coefficients of the CLUSTER, 
VOID, and FIELD samples are 0.98, 0.97, and 
0.96, respectively. The parameters of the linear 
regressions in the logarithmic plane are given in 
column (2) of Table 4. 
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Fig. 12. — Maximum circular velocity (a), and 
3D velocity dispersions averaged within Rh (b)and 
O.lRh (c) vs. halo mass for the Claii, Clf, VOID 
(Vail + Vi'ff), and FIELD samples. Sohd and 
dashed lines are linear regressions to the CLUS- 
TER and VOID samples, respectively. 



Our results confirm the tightness of the V^max- 
Mh relation, in this case, for different environ- 
ments. The relation is steeper for the VOID 
and FIELD samples than for the CLUSTER 
one. This is in agreement with the fact that the 
concentration-mass dependences for the former 
are shallower (almost absent) than for the lat- 
ter (see below). In Avila-Reese et al. (1999) 
it was shown how Vmax does depend on Mh; 
for example, for the NEW model (their eq. 8) 

Knax oc Afvir^^^5(cNFw), whcrc ^(cnfw) IS a func- 
tion of Cnfw and Cnfw depends on Mvir. In the 
limiting case of no dependence of Cnfw on Mvir, 
Kiiax oc A/vir^^"^, which is close to our results for 
VOID and FIELD halos. The steeper the depen- 
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Table 4 

Linear regression parameters (Kiax and rms velocities) 



Sample IgKiax = a + 61gMfc Igaj = c + cflgM^ IgtJc = e + /IgM^ 

aba c d a e f a 

z = 

CLUSTER -1.25 0.297 0.05 -1.52 0.320 0.08 -1.39 0.312 0.06 

VOID -1.76 0.336 0.03 -1.67 0.326 0.03 -1.60 0.328 0.04 

FIELD -1.76 0.334 0.03 -1.71 0.327 0.04 -1.54 0.322 0.04 

z = 1 

CLUSTER -1.55 0.324 0.04 -1.74 0.343 0.07 -1.43 0.319 0.06 

VOID -1.53 0.321 0.02 -1.52 0.321 0.02 -1.48 0.323 0.04 

FIELD -1.77 0.340 0.03 -1.82 0.370 0.05 -1.32 0.311 0.05 



dence of cnfw on Mvir, the shallower the 14iax- 
Mvir relation. 

As seen in Fig. 12(a), for a given mass, halos 
in voids and in the field have smaller Knax than 
halos in clusters. At IO^/i-IMq (3 IQi^/j-iMo), 
the circular velocity of VOID (FIELD) halos is 
typically ~ 20% (~ 12%) smaller than the cor- 
responding to CLUSTER halos. On the other 
hand, the CLUSTER sample shows the largest 
scatter in the Vmax-Mh relation. Several authors 
who modeled disk galaxy formation and evolution 
have shown that the TuUy-Fisher relation (TFR) 
of disk galaxies is well explained mainly as an im- 
print of the CDM halo V^ax-Mh relation (e.g., 
Mo et al. 1998; Steinmetz & Navarro 1999; Fir- 
mani & Avila-Reese 2000; Avila- Reese & Firmani 
2000; Buchalter et al. 2001). Thus, our results 
would imply that the TFR in different environ- 
ments could be slightly different (see also §6). 

In Fig. 13(a) we present the plot correspond- 
ing to Fig. 12(a) but at z = 1. For comparison 
purposes, the same linear regresions shown in Fig. 
12(a) (solid and dashed thin lines) are plotted. 
The V-axax^Mh relations in the different environ- 
ments are more similar at z — I than at z = 0. 
For a given mass, halos at z = 1 have on average 
slightly larger V^ax than at z = 0, specially for the 
VOID halos. This is because in a hierarchical sce- 
nario M/j increases with time typically more than 
Vmax does, shifting the relation to the high mass 
(low velocity) side. Nevertheless, the hierarchical 
mass growing can be stopped or even reversed in 
the cluster environment. The parameters of the 



linear regresions for the CLUSTER, VOID, and 
FIELD samples at z = 1 are given in column (2) 
of Table 4. The scatter of the Vmax-Mh relation is 
smaller at z = 1 than at z = 0, specially for the 
CLUSTER sample. 

5.2. Velocity dispersion vs mass 

Panels (b) and (c) of Figs. 12 and 13 show 

mass-weighted halo velocity dispersions vs mass 
for all samples at z = and z = 1, respectively. 
The total and central velocity dispersions are de- 
fined as follows: 

2 _ 47r/o^'* a^{r)p{rydr 
~ M(< Rh) 

2 _ 47r/° '^'' (J^{r)p{r)r''dr 
''^ ~ M{<Q.lRu) ' ^ ' 

where a{r) is the 3-dimensional velocity disper- 
sion at the spherical shell with radius r. The lin- 
ear fitting paranietc;rs for all the samples are given 
in third and fourth columns of Table 4. The be- 
haviour of the velocity dispersion-mass relations at 
z = for all the samples (environments) is similar 
to the one of the maximum circular velocity-mass 
relation (see also Antonucci-Delogu et al. 2001), 
reflecting this that most of the halos are close to 
their equilibrium state. In the plots for at vs M/j 
there are several outliers in the CLUSTER sam- 
ple. These are mostly truncated halos. They are 
not outliers in the ac vs M/, plot. Halos in voids 
and in the field have systematically smaller veloc- 
ity dispersions than halos in clusters, although the 



19 



1000 




Mass (h-iMg) 

Fig. 13. — Same as in Fig. 11 but for a snapshot 
at z — 1. For comparative purposes, the linear re- 
gressions of the CLUSTER and VOID samples at 
z = (Fig. 11) are shown (thin solid and dashed 
lines, respectively). 

difference becomes very small when the central ve- 
locity dispersion is used. 

6. The halo-galaxy connection 

Our study shows that present-day halos in the 
cluster environment are on average more con- 
centrated, more spherical, disaligned and rotate 
slower than halos in the void or field environ- 
ments. Our study also suggests differences in the 
evolution history of halos depending on the en- 
vironment: halos in dense environments assemble 
their masses apparently earlier than halos in low- 
density regions (see also Sheth & Tormen 2004). 
The question is whether all these differences in the 
halos produce visible differences in the luminous 



galaxies formed inside them. 

A common method in galaxy modeling is to 
build-up a luminous galaxy linked directly to a 
present-day CDM halo. Several numerical and 
analytical works show indeed that the present- 
day properties of the halos keep some memory of 
their evolution history (typical formation epoch, 
average shape of the MAH, major mergers, etc.). 
For example, the z = halo concentration de- 
pends on the shape of the MAH (Avila-Reese et 
al. 1998; Firmani & Avila-Reese 2000; Wechsler 
et al. 2002; Zhao et al. 2003), in such a way 
that fixing the concentration at z = 0, the main 
feature of the halo MAH (and of the correspond- 
ing baryon matter infall history) remains roughly 
determined. However, a direct {z — halo)- 
galaxy connection may fail for subhalos (specially 
for the massive ones and in cluster regions). The 
tumultous history of galaxy-size subhalos, since 
they infall in the cluster halo, affects in a com- 
plex way their present-day properties and distri- 
butions. The main effect is that of tidal stripping, 
or even total disruption, of subhalos as they pen- 
etrate into the cluster potential (e.g., Ghigna et 
al. 1998, 2000; Colfn et al. 1999, 2000; Taffoni et 
al. 2003; De Lucia et al. 2004; Kravtsov et al. 
2004a,b; Diemand et al. 2004; Gao et al. 2004a; 
Reed et al. 2005). Thus, the halo mass func- 
tion and the halo-to-galaxy mass ratio of galax- 
ies in clusters is expected to change dramatically 
with time. The study of galaxy evolution in clus- 
ters, at least in their inner regions, requires a full 
treatment of halo evolution and baryonic physics 
(Springel et al. 2001; Diaferio et al. 2001; Gao et 
al. 2004b). 

Following, we will carry out a very preliminar 
exploration of the effects on the disk galaxy prop- 
erties when varying the present-day CDM halo 
concentration and A, according to the results ob- 
tained above for different environments. For this, 
seminumerical models of disk galaxy evolution 
(Avila-Reese & Firmani 2000; Firmani & Avila- 
Reese 2000) will be used (see also van den Bosch 
2000). These models include in a self-consistent 
way the processes of formation and evolution of 
a spherical CDM halo and a disk in centrifugal 
equilibrium inside it, the adiabatic gravitational 
contraction of the halo due to disk formation, self- 
regulated SF and feedback, secular bulge forma- 
tion and other evolutionary processes. Most of the 
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Table 5 

Main properties of simulated disk galaxies 





Env.'' 


cnfw 


A 






B-V 


b/d 




Ms 




IO^Mq 








Mq/pc2 






10-2 


kpc 


IO^Mq 


km/s 


1 


V 


12 


0.048 


1.75 


0.64 


0.54 


1.1 


1.1 


3.61 


55 




CI 


19 


0.035 


2.42 


0.39 


0.74 


12.3 


0.7 


6.18 


71 


100 


V 


10 


0.048 


2.30 


0.50 


0.53 


8.8 


6.3 


512 


212 




CI 


13 


0.035 


2.83 


0.27 


0.72 


24.5 


3.8 


730 


278 



'^Environment: V=void, Cl=cluster 



z = Q galaxy properties depend mainly on (i) the 
present-day Mvir, (ii) the halo MAH (that deter- 
mines the concentration), (iii) the spin parameter 
A, and (iv) the disk mass fraction fd (= Md/M^i^). 
A brief description of the main ingredients of the 
model is presented in the Appendix. We notice 
that the A parameter in this case follows its stan- 
dard definition (e.g., Peebles 1969). This A pa- 
rameter is larger than the A' one measured for 
halos in §4 by factors typically of 1.25-1.50, de- 
pending on the halo concentration (Bullock et al. 
2001b). 

Wc model disk galaxies of a given baryon mass 
(Md) in halos with different concentration and 
spin parameters, emulating void and cluster en- 
vironments. The scminumcrical method is well 
suited for isolated (field and void) galaxies. For 
"void" galaxies, we assign the (high) median A 
(calculated from A') found for VOID halos and 
fix the z = Mvir {= Md/fd)- We then select, 
from random realizations, MAHs for this mass so 
that they yield the typical (low) concentrations 
of VOID halos at z = 0; these MAHs are ex- 
tended, implying late halo assembling. We care- 
fully choose roughly regular MAHs, without dra- 
matic changes in their shapes. The emulation of 
"cluster" galaxies with our method is more diffi- 
cult. We will assume that the (sub)halo was ac- 
creted into the cluster halo at z = 0.3 (3.4 Gyr ago, 
for the cosmology used here). Ours and previous 
simulations show that most subhalos in present- 
day cluster halos were accreted recently (~ 70% 
after z = 0.5), and that since then, these halos 
have lost typically 30-50% of their masses due to 
tidal striping (e.g., De Lucia et al. 2004; Nagai & 



Kravtsov 2005; Zentner et al. 2005; van den Bosch 
et al. 2005). Thus, we fix Mvir at 2; = 0.3 rather 
than at 2 = 0, and we assign a low A, typical of 
our CLUSTER halos at z = 0.3, which is slightly 
larger than at z = 0. Regarding concentrations, 
we fix them to the typical (high) concentrations of 
our CLUSTER halos by selecting the appropiate 
MAHs. These MAHs and the fact that the halo is 
fixed at 2; = 0.3, imply an early assembling of the 
halo/disk system. Since z < 0.3 the disk does not 
accretes more gas but it continues evolving. 

Models for two disk (baryonic) masses are cal- 
culated, Md = IO^Mq and IO^Mq. For each 
mass, we calculate two galaxy models correspond- 
ing roughly to the extreme environments, void and 
cluster regions, as explained above. We set the 
disk mass fraction in all cases to be /d = 0.03. The 
cosmologycal parameters are the same ones used 
in our N-body simulations (§2). Table 5 summa- 
rizes the main halo and disk input parameters as 
well as the obtained properties of our simulated 
galaxies. Recall that in the case of the "cluster" 
galaxies, the halo growth and gas infall are trun- 
cated at z = 0.3 to account for the fact that it 
was accreted into the cluster at this time. The 
late-accreted halos are typically located in the pe- 
riphery of the cluster halo (e.g., De Lucia et al. 
2004; Nagai & Kravtsov 2005) and are expected 
to host mostly spiral or SO galaxies. 

Looking at Table 5, one sees that the depen- 
dence of halo concentration and A on environment 
found in this work should produce some changes 
in the disk galaxy properties. As previously re- 
ported (Avila- Reese & Firmani 2000), A influ- 
ences mainly the disk scale length and surface 
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brightness, the gas fraction and the secular bulge- 
to-disk ratio, while the MAH fixes the halo con- 
centration and influences mainly the SF history, 
galaxy integral color and the scatter of the TFR. 
Thus, our expectation is that disk galaxies of a 
given mass formed in low-A and highly con- 
centrated halos, with a gas infall history truncated 
early (cluster environment), are preferentially of 
earlier morphologycal types, redder, shifted to the 
higher velocity side in the TFR, and have higher 
surface brightness, smaller scale lengths, and lower 
gas fractions than disk galaxies formed in high-A 
and low-concentration halos (void environment). 
All these trends are namely seen in Table 5. The 
stellar central surface density, Sq, and disk scale 
radius, r^, are the parameters of the exponential- 
law fit to the model stellar disk. In fact, the disks 
formed in CDM halos typically are not exactly ex- 
ponentials, but their surface density profiles are 
more concentrated in the center and with an ex- 
cess in the periphery compared with the exponen- 
tial law (Firmani & Avila- Reese 2000; Bullock et 
al. 2001b). 

Our results show that owing to only differ- 
ences in the halo properties (related to the en- 
vironment), the galaxy disks can have significant 
differences. According to Table 5, disks formed 
in halos with properties typical of halos in voids 
are ~ 0.2 mag bluer in {B — V) color, have cen- 
tral stellar surface densities lower by 4 times 
(1.5 mag/arcsec^), and gas fractions higher by 1.5- 
2 times than disks formed in halos with proper- 
ties typical of the cluster-periphery environment. 
According to our models, the low mass disks 
have lower central surface brightnesses and higher 
gas fractions than the high mass ones, while the 
{B — V) colors do not change significantly with 
mass. The (secular) bulge-to-disk ratio depends 
strongly on mass, being this ratio larger for mas- 
sive disks. There is also a significant increasing 
of this ratio from halos typical of voids to those 
typical of the cluster periphery. This implies that 
part of the morphology- density relation is due to 
changes in the properties of the galaxy halos with 
environment. 

As mentioned in §5.1, one expects also differ- 
ences with environment in the TFR. The slope 
of the model stellar TFR is ~ 3.4 (Firmani & 
Avila-Reese 2000). By using this slope, we cor- 
rect the velocities due to the variations in the 



stellar masses of each one of the models with 
Md = IO^Mq and Md = 10"Mq presented in 
Table 5. We obtain that at Md = IO^^Mq, the dif- 
ference in Vrnax for our void and cluster-periphery 
galaxies is about 20%; the difference in mass, ex- 
pressed in magnitudes, is 0.6 mag, i.e. we predict 
that the zero-point of the TFR of void galaxies 
at Md = 10"Mq should be brighter by ~ 0.6 
mag than the one of cluster-periphery galaxies. At 
Md = IO^Mq, the predicted difference is ~ 0.4 
mag. 

Real galaxies show likely more pronounced dif- 
ferences in their observational properties as a func- 
tion of the environment than those obtained with 
our models (see §1.1 for references). Thus, other 
physical processes not considered here should cer- 
tainly play an important role in galaxy dynam- 
ics and evolution. For example, we did not treat 
in detail the angular momentum distribution mis- 
alignment and did not take into account the halo 
triaxiality and the shape-to-rotation axis misalign- 
ment. These properties change with environment 
as was shown in §§3 and 4. Besides, deep in the 
cluster regions, the external effects such as tidal 
and ram pressure stripping, strong interactions at 
early epochs, and galaxy harassment (Moore et 
al. 1996), are likely the dominant ones in shaping 
morphology and other galaxy properties (see the 
references above). 

Finally, we note that most of the differences in 
the halo properties with environment seen at z = 
dissapear at z = 1, although the scatters in the lat- 
ter epoch are very large. As mentioned in §1.1, ob- 
servations also show that the morphology-density 
relations flattens at higher redshifts. The main 
morphological evolution is seen for the galaxies in 
the high-density environment. Halos in clusters 
are also those that change most their properties. 

7. Summary and Conclusions 

Observations show that morphology and sev- 
eral properties of galaxies change systematically 
with environment. According to the current 
paradigm of galaxy formation, galaxies assem- 
ble inside CDM halos. We state then two ques- 
tions: (i) do the properties of galaxy-size CDM 
halos/subhalos change systematically with envi- 
ronment?, and if this the case, (ii) do these changes 
in the halo/subhalo properties affect the luminous 
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galaxies in the direction that observations show? 

We studied the first question by means of 
high-resolution ACDM cosmological N-body sim- 
ulations. Several clusters with their surround- 
ings (CLUSTER sample) and void regions (VOID 
sample) were selected from large-box simulations, 
and they were resimulated with high resolution. 
We also use other simulations, in particular two 
60/i~^Mpc-box simulations from which the FIELD 
sample of parent halos was extracted. The local 
density contrast in our CLUSTER, VOID and 
FIELD samples are on average larger than 34, 
smaller than -0.8 and roughly 0, respectively. The 
second question, in a first attempt, was disscused 
by using standard seminumerical models of disk 
galaxy evolution. Disk galaxies were modeled 
inside CDM halos with properties that we have 
found they have in the different environments. 
Following, we summarize our main results: 

• For masses < 5 x lO^^h^^MQ, halos from 
the CLUSTER sample are on average ~ 40% 
more concentrated and have ^ 2 times higher cen- 
tral densities p_2 than halos in voids at z = 0. 
While for halos in cluster regions the concentra- 
tion parameters cnfw, Ci/5 and c^, and the density 
P-2 decrease on average with mass, for halos in 
voids these concentrations and p_2 do not seem to 
change with mass. The slope of the cnfw — -^vir, 
Ci/5 — Mvir and c^ — Mh dependences for the for- 
mer sample is ^ —0.1. In the mass range of the 
parent FIELD halos analyzed here, concentrations 
and /3_2 are also smaller on average than those of 
halos in clusters. The concentrations and p^2 of 
FIELD halos decrease on average with mass but 
less rapid than halos from the CLUSTER sample. 
The scatters of all of these parameters are larger 
for CLUSTER halos than for halos in less dense 
environments. For example, for reasonably well 
fitted NFW halos from the CLUSTER and VOID 
samples, A(logCNFw)~ 0.13 and w 0.09, respec- 
tively. 

• All the differences mentioned above become 
less pronounced when comparing only the par- 
ent halos from the different samples (CLUSTER: 
~ 60%; VOID: ~ 95%; FIELD: all). Therefore, 
the CLUSTER halos are more concentrated and 
internally denser than the VOID and FIELD ones 
due partially to a local halo-subhalo effect (sub- 
halos are a significant fraction of the CLUSTER 
sample, and subhalos are systematically more con- 



centrated than their parent halos, see also, e.g., 
Bullock et al. 2001a). However, we find that the 
parent halos and subhalos from the CLUSTER 
sample are still significanlty different on average 
than the parent halos and subhalos from the VOID 
and FIELD samples. Therefore, the differences in 
halo properties are certainly also due to a pure 
global environmental effect, related mainly to the 
typical halo formation epoch: halos in dense envi- 
ronments assemble their masses earlier than halos 
in low- density regions. 

• Halos in dense environments are more spheri- 
cal than halos in less dense environments at z = 0. 
The minor- to-major axis ratios of CLUSTER ha- 
los are on average ~ 1.2 times lower than those of 
the FIELD halos. For a given epoch, the elliptic- 
ity of the halos tends to increase with mass. For 
a given mass, the ellipticity changes with age, the 
younger halos having on average larger ellipticities 
than the older halos. 

• The spin parameter of CLUSTER halos is on 
average 1 .3 — 1 .4 times lower than the one of VOID 
or FIELD halos at z = 0. This is likely a con- 
sequence of both global (environmental) and lo- 
cal effects, i.e. the dependence of halo formation 
epoch on environment, and the tidal stripping and 
"tumultuos" histories that halos suffer in locally 
high-density regions, respectively. We do not find 
significant differences in the A' distribution of par- 
ent halos and subhalos from the CLUSTER sam- 
ple. The CLUSTER halos appear to have a less 
aligned intrinsic angular momentum distribution 
than the VOID and FIELD ones, but this result 
needs to be confirmed by future analysis because 
of the large errors involved in the determination 
of cos^^x/2- The spin parameter docs not change 
significantly with mass in any environment and its 
distribution is well approximated by a lognormal 
function with a larger width and a lower peak A' in 
cluster regions than in voids and the field. The an- 
gular momentum axis of halos tends to be aligned 
with their minor principal axis, this behaviour be- 
ing more common in the less dense environments 
than in the high-density ones. 

• A tight V^nax — Mh relation is seen for halos 
in all the environments. The relation is shallower 
(slope of 0.30) and more scattered for clusters and 
their surroundings than for the void and field envi- 
ronments (slopes of 0.33 — 0.34). If this relation 
is the basis of the observed TFR, then slight dif- 
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ferences in the TFR are expected in different envi- 
ronments. The Vmax — Mh relation sUghtly shifts 
to the lower mass side at ^; = 1 in all the envi- 
ronments, but in particular in the voids. Similar 
results to the Knax — relation are obtained for 
the (Tc— and at — Mh relations, showing that most 
of the halos are close to their equilibrium state. 

• The differences in halo properties with en- 
vironment seen at 2; = drastically diminish at 
z = 1. Interestingly enough, a similar result was 
found for the observed morphology-environment 
relation. As expected, the concentration parame- 
ters decrease as we go from z = to z = 1, being 
this change more pronounced for the CLUSTER 
halos, particularly for the less massive ones. The 
spin parameter of CLUSTER halos at z = 1 is on 
average significantly larger than at z = 0, while 
for less dense environments, the spin parameter at 
z = 1 is on average only slightly larger than at 
z = 0; thus, A' evolves in a different way for halos 
in cluster-like regions than for halos in less dense 
environments. Halos in the void and field envi- 
ronments are systematically less aligned at z = 1 
than at z = 0, while halos in the cluster regions do 
not show any systematical change with redshift in 
a statistical sense. 

• Disk galaxies modeled in a self-consistent 
fashion inside ACDM halos with the present-day 
concentrations and spin parameters found here for 
halos in the different environments present sys- 
tematical differences: the galaxies formed in ha- 
los typical of cluster-periphery environment have 
higher surface density, circular velocity and sec- 
ular bulge-to-disk ratio, lower gas fraction, and 
are redder than for those formed in halos typical 
of void environment. These trends agree quali- 
tatively with observations but are not enough to 
explain the observed differences of galaxy proper- 
ties with environment. We predict that the TFR 
of galaxies in low and high-density environment is 
different. 

From our study we conclude that most of the 
properties of galaxy-size halos at z = change 
with environment in a statistical sense, the largest 
differences in the structural properties being for 
the less massive halos (sub-L* galaxies). The 
main changes with environment occured after 
z ~ 1 and the most affected halos are the sub- 
halos in the CLUSTER sample. The differences 
in halo concentrations and spin parameters along 



the different environments influence on the prop- 
erties of galaxy disks formed inside these halos 
and in the same direction that observations show. 
However, the inclusion of angular momentum 
dissalignment, triaxiality, and shape-rotation axis 
dissalingment -halo properties that also change 
with environment - in the models of disk galaxy 
evolution is necessary in order to attempt to re- 
produce in more detail the observed relations of 
disk galaxy properties with environment. The 
halo merging history and its influence on lumi- 
nous galaxies should also be taken into account, 
in particular for modelling galax;ies in the high- 
density environments. 

The halo properties discussed above and their 
changes with environmnet are ultimately related 
to the initial cosmological conditions. Based on 
our results, we claim that the observed galaxy 
properties-environment relations are partially es- 
tablished by the initial cosmological conditions 
(nature), in the sense that galaxy halos in more 
dense environments assemble earlier. However, 
our results point out that the external astrophys- 
ical factors (nurture) should play also an impor- 
tant role in the observed trends of changing galaxy 
properties with environment, specially in the high- 
est density regions (clusters). 

Computer simulations presented in this paper 
were done at the Leibnizrechenzentrum (LRZ) 
in Munich and at the John von Neumann Insti- 
tute for Computing Jiilich. We acknowledge the 
anonymous referee whose helpful comments and 
suggestions improved several aspects of the pa- 
per. We are also grateful to Ricardo Florcs for 
kindly providing a copy of his program to com- 
pute the halo ellipticities. This work has been sup- 
ported by a bilateral CONACyT-DFG (Mexico- 
Germany) grant, and by CONACyT grants 36584- 
E and 40096-F. 



APPENDIX 

Here we present the main physical ingredients 
of the self-consistent evolutionary models used in 
§6. For details sec Firmani & Avila-Recsc (2000); 
Avila-Reese & Firmani (2000). The disk is built 
up within a growing ACDM halo. An extended 
Press-Schechter approach (Lacey & Cole 1993; 
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Avila-Rccsc ct al. 1998) is used to generate the; 
statistical MAHs of the halos from the primor- 
dial density fluctuation field, and a generalized 
secondary infall model is applied to calculate the 
time by time virialization of the accreting mass 
shells (Avila- Reese et al. 1998). The evolution 
and structure of the ACDM halos calculated this 
way agree well with results from cosmological N- 
body simulations (Avila- Reese et al. 1999). Halos 
assembled through early active MAHs end more 
concentrated on average than halos with extended 
MAHs. 

The halo mass shells are assumed to have 
aligned rotation axis with specific angular momen- 
tum given by jsh{tv) = dJ{ty)/dMy{ty), where 
J = XGMy^'^/ |£;|^/^ J, M„ and E are the total 
angular momentum, mass and energy of the halo 
at the shell virialization time ty. The spin param- 
eter, A, is assumed to be constant in time. As 
the result of the assembling of these mass shells, a 
present day halo ends with an angular momentum 
distribution close to the (universal) distribution 
measured by Bullock et al. (2001b) in N-body 
simulations. A fraction fd of the mass of each shell 
is assumed to cool down and form a disk layer in 
a dynamical time. The radial mass distribution 
of the layer is calculated by equating its specific 
angular momentum to that of its final circular 
orbit in centrifugal equilibrium (detailed angular 
momentum conservation is assumed). The super- 
position of these layers form the disk. 

The gravitational interaction of disk and inner 
halo is calculated using a generalized adiabatic in- 
variance formalism. This formalism differs from 
the usual one (e.g., Flores et al. 1993) in that 
we take into account the ellipticity of the orbits 
(the circular orbit assumption is relaxed); recent 
numerical studies (Gnedin et al. 2004) confirm 
the validity of our approach. 

The disk SF at a given radius (azhimutal sym- 
metry is assiuncid) is triggered by the Toomre 
gas gravitational instability criterion and self- 
regulated by a vertical disk balance between the 
energy input due to SNc and the turbulent energy 
dissipation in the ISM. This physical prescription 
for disk SF yields naturally a Schmidt law with an 
index n < 2, slightly varying along the disk (Fir- 
mani Hernandez & Gallagher). The SF efficiency 
depends on the gas surface density determined 
mainly by A, and on the gas accretion rate de- 



termined by the cosmological MAH. Finally, we 
estimate the mass of the (secularly formed) bulge 
as the inner disk mass where the Toomre stellar 
parameter indicates disk instability (see also van 
den Bosch 1998) 
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